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Summary. Axion emission by hot and dense plasmas is a new energy-loss channel 
for stars. Observational consequences include a modification of the solar sound-speed 
profile, an increase of the solar neutrino flux, a reduction of the helium-burning 
lifetime of globular-cluster stars, accelerated white-dwarf cooling, and a reduction 
of the supernova SN 1987A neutrino burst duration. We review and update these 
arguments and summarize the resulting axion constraints. 



1 Introduction 

The "outer space" of astrophysics and cosmology provides a natural labora- 
tory for the "inner space" of elementary particle physics. Usually one may first 
think of the early universe or perhaps high-energy cosmic rays for arguments 
in favor or against a new particle-physics model. However, in the case of axions 
the low energies available in stars are well suited for very sensitive tests. 

The basic idea is very simple. Stars are powerful sources for weakly inter- 
acting particles such as neutrinos, gravitons, hypothetical axions, and other 
new particles that can be produced by nuclear reactions or by thermal pro- 
cesses in the stellar interior. Even when this particle flux cannot be directly 
measured, the properties of stars themselves would change if they lost too 
much energy into a new channel. This "energy-loss argument" has been widely 
used to constrain a long list of particle properties, see [1-7] for early exam- 
ples and [8-12] for extensive reviews. We summarize here the main arguments 
that have been put forward, the observational evidence, and the resulting 
constraints for "invisible axions." 

To this end we review in Sect. [2] the axion interactions with photons and 
fermions. In Sect. [3] we consider the Sun as an axion source, notably by the 
Primakoff process, and review limits on the axion-photon interaction strength 
by helioseismology, the measured neutrino flux, and the CAST experiment. 
In Sects. 131 El and [HI we review axion limits from globular cluster stars, white- 
dwarf coohng, and supernova SN 1987A, respectively. We summarize these 
constraints in Sect. [7] in juxtaposition with cosmological arguments. 
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2 Axion Interactions 



The particle-physics motivation for "invisible" axions and their main proper- 
ties were introduced in Chap. 1 of this volume. Before turning to their role 
in stars, we briefly review the phenomenological properties of these pseudo 
Nambu-Goldstone bosons of the Peccei Quinn (PQ) symmetry. The mass and 
interaction strength with ordinary particles is approximately given in terms 
of the relevant 7r° properties, scaled with J-^jja where = 92MeV is the 
pion decay constant and /q is the PQ scale or axion decay constant. The 
normalization of /a is defined by the axion-gluon interaction 
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where a is the axion field, G the gluon field-strength tensor, G its dual, and b 
a color index. Color anomaly factors have been absorbed in this definition of 
fa, which is the quantity that is relevant for all low-energy phenomena [13]. 

The PQ symmetry is explicitly broken at low energies and axions acquire 
a small mass. Unless there are non-QCD contributions, perhaps from Planck- 
scale physics [14, 15], the mass is 

_ Um^ _ 6.0 eV 

""""l + z /. -/JlOeCeV 

where z = niu/rnd is the up/down quark mass ratio. We will follow the previ- 
ous axion literature and usually assume the canonical value z — 0.56 [16, 17], 
although it could vary in the range z = 0.3-0.6 [12]. 

Another generic property of axions is their two-photon interaction that 
plays a key role for most searches, 

-^177 — i-^L F^yF'^'^a = —gajf E • B a . (3) 

Here, F is the electromagnetic field-strength tensor, F its dual, and E and B 
the electric and magnetic fields, respectively. The coupling constant is 

_ a f E 24 + z\_al'E 2 4 + z\ 1 + z nia 
^"^"'"2^^iV"3TT^j " 2^^]v"3TT^j7i7^^^' ^ ^ 

where E and N, respectively, are the electromagnetic and color anomaly of the 
axial current associated with the axion field. E/N = 8/3 in grand unified mod- 
els, e.g. the DFSZ model [20,21], whereas E/N = in the KSVZ model [18,19]. 
While these cases are often used as generic examples, in general E/N is not 
known so that for fixed fa a broad range of gaj-y values is possible [22]. Still, 
barring fine-tuned cancelations, 5077 scales from the corresponding pion in- 
teraction by virtue of the relation ([4]). Taking the model-dependent factors to 
be of order unity, this relation defines the "axion line" in the ma~ga-yy plane. 
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Axions or axion-like particles with a two-photon vertex decay into two 
photons with a rate 
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where the first expression is for general pseudoscalars, the second applies 
specifically to axions, and the numerical one assumes z = 0.56 and the 
hadronic case E/N — 0. Comparing with the age of the universe of 4.3 x 10"'^'' s 
reveals that axions decay on a cosmic time scale if iria > 20 eV. 

The interaction with fermions j has a derivative structure so that it is 
invariant under a — > a + ao as behooves a Nambu-Goldstone boson, 

C C 'TTt 

Lajj = ^^j7''75^jd^,a or - \ ^p-^j^^^^a . (6) 

^/a J a 

Here, is the fermion field, rrij its mass, and Cj a model-dependent numer- 
ical coefficient. The combination Qajj = Cj-mj/ fa plays the role of a Yukawa 
coupling and aajj = g'^jj/^'K of a "fine-structure constant." The pseudoscalar 
form is usually equivalent to the derivative structure, but one has to be careful 
in processes where two Nambu-Goldstone bosons are attached to one fermion 
line, for example an axion and a pion attached to a nucleon in the context of 
axion emission by nucleon bremsstrahlung [23,24]. 

In hadronic models such as KSVZ [18,19], axions do not couple to ordinary 
quarks and leptons at tree level, whereas in the DFSZ model [18, 19] 

C^e = ^. (7) 

Here, cot (3 is the ratio of two Higgs vacuum expectation values of this model. 

For micleons, the dimensionless couplings Cn,p are related by generalized 
Goldbcrgcr-Treiman relations to nucleon axial- vector current matrix elements, 

Cp = (C„ -r])Au+ {Cd - r]z) Ad + [C, - rjw) As , 

Cn = {Cu -ii)Ad+ {Cd - iiz) Au + {Cs - vw) As . (8) 

Here, rj = {1 + z + w)~^ with z = niu/md and w = mu/rris <?C z. The 
quantities Aq represent the axial-vector current couplings to the proton, 
AqS^ = (p|^7p75(?|p} where is the proton spin. 

Neutron beta decay and strong isospin symmetry tell us that Au — Ad = 
F + D = 1.267 ±0.0035 whereas hyperon decays and flavour SU(3) symmetry 
imply Au + Ad - 2As = 3F - D = 0.585 ± 0.025. Recent determinations of 
the strange-quark contribution are As = —0.08 ± O.Olgtat ± 0.05syst from the 
COMPASS experiment [25] and As = -0.085±0.008exp±0.013theor±0.009evoi 



4 Georg G. Raffelt 



from HERMES [26] , in agreement with each other and with an early estimate 
of As = -0.11 ± 0.03 [27]. We thus adopt the estimates 

Au = +0.841 ±0.020 , 
Ad^ -0.426 ± 0.020 , 

As = -0.085 ±0.015 , (9) 

that are very similar to previous values used in the axion literature. 

The uncertainty of the axion-nucleon couplings is dominated by the large 
uncertainty of z = 0.3-0.6 that was mentioned above. For hadronic axions 
we have Cu,d.s = so that Cp — —0.55 and C„ = +0.14 for z = 0.3 and 
Cp — —0.37 and C„ = —0.05 for z — 0.6. Therefore, while it is well possible 
that Cn — 0, Cp does not vanish within the plausible z range. In the DFSZ 
model we have C„ = | sin^ (i and Cd = \ cos^ (3. Even with the large allowed z 
range, Cn and Cp never vanish simultaneously. An extreme case is cos^ /3 = 
where Cp = for z = 0.3, but in this case Cn — —0.27. 



3 The Sun as an Axion Source 

3.1 Axion Flux from the Primakoff Process 

The Sun would be a powerful axion source. This flux can be searched directly, 
notably by the CAST experiment. Its sensitivity is competitive with the glob- 
ular cluster limits (Sect. |4|) for hadronic models. In this case the dominant 
emission process is the Primakoff effect [28], i.e., particles with a two-photon 
vertex transform into photons in external electric or magnetic fields. There- 
fore, stars produce axions from thermal photons in the fluctuating electro- 
magnetic fields of the stellar plasma [6] . 

Calculating the solar axion flux is straightforward except for the proper 
inclusion of screening effects [29,30]. The transition rate for a photon of energy 
E into an axion of the same energy (recoil effects are neglected) is [31] 

(-J.)-(-f)-]' 

where T is the temperature (natural units with h = c = = 1 are used). 
The screening scale in the Debye-Hiickel approximation is 

nuclei 

where rie is the electron density and rij that of the j-th ion of charge Zj. 
Near the solar center Kg ~ 9keV. Note that (ks/T)^ k, 12 is nearly constant 
throughout the Sun whereas it is about 2.5 throughout the core of a low-mass 
helium-burning star. 



~ 32^ 
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Ignoring the plasma frequency for the initial-state photons, the energy-loss 
rate per unit volume is [10, 31] 

Q^-^F^ (12) 

where F is a numerical factor of order unity. For (ks/T)^ = 2.5 and 12 one 
finds F = 0.98 and 1.84, respectively. 

Integrating over a standard solar model, one finds an axion flux at Earth 
that is is well approximated by {E in keV) 

^ = gl^ 6.0 X IQi" cm-2 s"! keV-^ iJ^ ^si g-^^/i-^os _ (13) 

where qiq — 5077/(10^^" GeV~^). The integrated flux parameters are 

= gl^ 3.75 X 10" cm~2 s"i , 
La = glo 1-85 X IO-^Lq . (14) 

The maximum of the distribution is at 3.0 keV, the average energy is 4.2 keV. 

3.2 Solar Age 

The properties of the Sun itself constrain this flux. The axion losses lead 
to an enhanced consumption of nuclear fuel. The standard Sun is halfway 
through its hydrogen-burning phase so that the solar axion luminosity should 
not exceed its photon luminosity Lq. 

As an example we recall that a magnetically induced vacuum dichroism 
observed by the P VLAS experiment [32] , if interpreted in terms of an axion- 
like particle (ALP), requires [32,33] 5^77 = 2-5 x lO^^GeV^^ and rria — 
1-1.5 meV. With this coupling strength, the Sun's ALP luminosity would 
exceed Lq by a factor of a million and thus could live only for about 1000 years. 
Perhaps this problem can be circumvented, but it is noteworthy that even a 
crude astrophysical argument severely constrains the particle interpretation 
of the PVLAS signature. 

3.3 Helioseismology 

For a more refined constraint we note that a model of the present-day Sun, 
with the integrated effect of axion losses taken into account, would differ from 
a standard solar model. The modified sound-speed profile can be diagnosed 
by helioseismology, providing a conservative limit [34] 

ffa77 < 1 X 10"^GeV"\ (15) 

corresponding to La < 0.20 L©. More recent determinations of the solar metal 
abundances have diminished the agreement between standard solar models 
and helioseismology [35], but these modifications do not change the limit 
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3.4 Solar Neutrino Flux 

The energy loss by solar axion emission requires enhanced nuclear burning 
and thus a somewhat increased temperature in the Sun. Self-consistent solar 
models with axion losses reveal that g^^^ = 4.5 X 10-10 GeV-1 implies a 20% 
increase of the solar neutrino flux [34]. For ga-y^ = 10 x 10-^" GeV-^ the 
increase would be a factor of 2.4. 

The measured all-flavor *B neutrino flux is 4.94 x lO^cm-^s""'^ with an 
uncertainty of about 8.8% [36,37]. The old standard solar model predictions 
were about 5.7-5.9 in the same units, whereas the new metal abundances 
imply 4.5-4.6, each time with a 16% "theoretical Icr error" [35]. Therefore, 
the measured neutrino fluxes imply a limit 

ffa77 < 5 X 10"i°GeV-\ (16) 

corresponding to La ^ 0.04 L©. A more precise limit with a realistic error 
budget would require self-consistent solar models on a finer spacing of ga-y-f 

3.5 Searches for Solar Axions 

The solar axion flux can be searched with the inverse Primakoff process where 
axions convert to photons in a macroscopic B field, the "axion helioscope" 
technique [38]. One would look at the Sun through a "magnetic telescope" 
and place an x-ray detector at the far end. The conversion can be coherent 
over a large propagation distance and is then pictured as a particle oscillation 
effect [39]. 

Early helioscope searches were performed in Brookhaven [40] and Tokyo 
[41,42]. Solar axions could also transform in electric crystal fields, but the lim- 
its obtained by SOLAX [43], COSME [44], and DAMA [45] are less restrictive 
and require a solar axion luminosity exceeding (|15p and (|16p . i.e., these limits 
are not self-consistent. 

The first helioscope that can actually reach the "axion line" is the CERN 
Axion Solar Telescope (CAST). The non-observation of a signal above back- 
ground leads to a constraint [46] 

ga77 < 1-16 X 10^1" GeV^i (95% CL) for < 0.02 eV . (17) 

For larger masses, the axion-photon transition is suppressed by the energy- 
momentum mismatch between particles of different mass. The full rate can be 
restored in a narrow range of masses by providing the photons with a refractive 
mass in the presence of a low-Z gas [47], a method that was already used in 
the Tokyo experiment [42] and is also used in the ongoing GAST Phase II. 
GAST is foreseen to reach eventually nia < 1 eV. 
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3.6 Do Axions Escape from the Sun? 

CAST can detect axions only if they actually escape from the Sun. Their 
mean free path (mfp) against the Primakoff process is the inverse of (fTO| . For 
4keV axions and with T k, 1.3 keV and Kg ~ 9keV at the solar center we find 
Ao » 9io & X 10^'* cm K. 8 x 10^^ i?©, or about 10~^ of the radius of the 
visible universe. Therefore, ga^-y would have to be more than 10^ larger than 
the CAST limit for axions to be re-absorbed in the Sun. 

Even in this extreme case they are not harmless because they would carry 
the bulk of the energy flux that otherwise is carried by photons. The mfp 
of low-mass particles in the trapping regime should be shorter than that of 
photons (about 10 cm near the solar center) to avoid a dramatic modification 
of the solar structure [48]. This requirement is so extreme that for anything 
similar to axions the possibility of re-absorption is not a serious possibility. 

4 Globular-Cluster Stars 

4.1 Helium-Burning Lifetime and the Axion-Photon Interaction 

A restrictive limit on 3077 arises from globular-cluster stars. A globular cluster 
is a gravitationally bound system of stars that formed at the same time and 
thus differ primarily in their mass. A globular cluster provides a homogeneous 
population of stars, allowing for detailed tests of stellar-evolution theory. The 
stars surviving since formation have masses somewhat below 1 Mq. In a color- 
magnitude diagram (Fig. [1]) , where one plots essentially the surface brightness 
vs. the surface temperature, stars appear in characteristic loci, allowing one 
to identify their state of evolution. 

The stars on the horizontal branch (HB) have reached helium burning, 
where their core (about 0.5 Mq) generates energy by fusing helium to car- 
bon and oxygen with a core-averaged energy release of about 80ergg~^ s^^. 
A typical density is lO'^gcm"^ and a typical temperature 10^ K. The Pri- 
makoff energy loss rate (|12p implies that the energy-loss rate per unit mass, 
e = Q/g, is proportional to T'^ / Q- Averaged over a typical HB-star core one 
finds ((T/IO^K)^ (lO^gcm^'^/g)) « 0.3. Therfore, the core-averaged energy 
loss rate is about g^Q 30ergg^^ s^^. The main effect would be accelerated 
consumption of helium and thus a reduction of the HB lifetime by a factor 
80/(80 + 30.g?o), i.e., by about 30% for = 1. 

The HB lifetime can be measured relative to the red-giant evolutionary 
time scale by comparing the number of HB stars with the number of RGB 
stars that are brighter than the HB. Number counts in 15 globular clusters [49] 
show that this number ratio agrees with expectations within 20-40% in any 
one cluster, where the error is mostly statistical because typically only about 
100 HB stars were present in the used fields of view. Compounding the results 
of all 15 clusters, the helium-burning lifetime agrees with expectations to 
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Fig. 1. Color-magnitude diagram for tlie globular cluster M3, based on 10,637 
stars [50]. Vertically is the brightness in the visual (V) band, horizontally the dif- 
ference between B (blue) and V brightness, i.e. a measure of the color and thus 
surface temperature, where blue (hot) stars lie toward the left. The classification 
for the evolutionary phases is as follows [51]. MS (main sequence): core hydrogen 
burning. BS (blue stragglers). TO (main-sequence turnoff): central hydrogen is ex- 
hausted. SGB (subgiant branch) : hydrogen burning in a thick shell. RGB (red-giant 
branch): hydrogen burning in a thin shell with a growing core until helium ignites. 
HB (horizontal branch): helium burning in the core and hydrogen burning in a shell. 
AGB (asymptotic giant branch): helium and hydrogen shell burning. P-AGB (post- 
asymptotic giant branch): final evolution from the AGB to the white-dwarf stage 



within about 10% [10,11]. Of course, with modern data these results likely 
could be improved. Either way, a reasonably conservative Umit is 

ga77 < 10~^"GeV"^ . (18) 

It is comparable to the CAST Hmit p?]) . but applies for higher masses. The 
relevant temperature is about 10 keV so that significant threshold effects begin 
only at about nia > 30 keV. For QCD axions the coupUng increases with mass 
so that the Hmit reaches to even larger masses. 

In the hehum-burning core, convection and semi-convection dredges he- 
lium to the burning site so that 25-30% of all helium is burnt during the HB 
phase. Therefore, while the standard theoretical predictions depend on a phe- 
nomenological treatment of convection, there is limited room for additional 
energy supply, even if the treatment of convection were grossly incorrect. 




Astrophysical Axion Bounds 9 



4.2 Helium Ignition and the Axion-Electron Interaction 

Stars on the red-giant branch (RGB) have a degenerate hchum core with a 
typical density 10^ gcm^"^ and T « 10'* K. Helium ignites at a critical combi- 
nation of Q and T. Therefore, helium ignition can be delayed by axion cooling. 
This implies that the core grows more massive before helium ignites. One con- 
sequence is that the RGB will extend to brighter stars, i.e., the brightness of 
the brightest red giant in a given globular cluster signifies the core mass at 
helium ignition. Detailed studies reveal that the core mass at helium igni- 
tion agrees with theoretical expectations within 5-10% [10,52-55]. In turn, 
this implies that a novel energy-loss rate at T = 10® K and an average den- 
sity {q) — 2 X lO^gcm"'^ should not exceed about 10ergg~^s~^. At these 
conditions the standard neutrino emission is about 4ergg~^ s~^. 

The helium-burning lifetime is useful to constrain the axion-photon inter- 
action because the Primakoff rate is suppressed in the degenerate red-giant 
cores and thus is more effective in HB stars. The helium-ignition argument, 
on the other hand, is useful when the emission rate is larger on the RGB than 
on the HB as for bremsstrahlung e -t- Ze — > Ze -\- e -\- a. For the conditions in 
a red-giant core one finds ebrcms ~ ctaee 2 x 10^'' ergg~^ s'^ [54] so that 

aaee < 0.5 X lO^^s or gaee < 3 X IQ-^^ . (19) 

In the DFSZ model this limit corresponds to fa/ cos^ /3 > 0.8 x 10^ GeV, 
ma < 9meV/cos^/3 and ga-y-ycos'^ (3 < 1.2 x 10~^^ GeV"^ 

4.3 Asymptotic Giant Branch (AGB) Evolution 

For axion-clcctron interactions near or even below the bound (|19p . the emis- 
sion will strongly affect the evolutionary behavior of AGB stars [56]. However, 
these results have not been linked closely enough to observational data to ob- 
tain new limits or discover evidence for axion emission. 



5 White-Dwarf Cooling 

The degenerate core of a low-mass red giant before helium ignition is essen- 
tially a helium white dwarf. After the HB phase, when helium burning has 
ended, low-mass stars once more ascend the red-giant branch as "asymptotic 
giants" (AGB stars) . They have a degenerate carbon-oxygen core and helium 
burning in a shell. Fast mass loss creates a "planetary nebula" surrounding 
a compact remnant, a white dwarf, that first cools by neutrino emission and 
later by surface photon emission. 

The observed white-dwarf luminosity function reveals that their cooling 
speed agrees with expectations, constraining new cooling agents such as axion 
emission [10,57-59]. The resulting limit on the axion-electron coupling of 
Oiaee ^ 1 X 10^^^ is Comparable to the globular-cluster limit of (fT9|) . 
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The cooling speed of individual white dwarfs can be estimated in some 
cases where they appear as ZZ Ceti stars, i.e., when they are pulsationally 
unstable and when the period decrease P /P can be measured, a quantity that 
is sensitive to the cooling speed. A well-studied case is the star G117-B15A. 
For some time it seemed to be cooling too fast, an effect that could have been 
attributed to axion cooling with aaee = 0.2-0.8 x 10~^^ [60]. More recent 
analyses no longer require a new cooling channel, allowing one to set a limit 
on axion losses corresponding [61,62] 

aaee < 1-3 X lO'^^ Or gaee < 1.3 X 10"" (20) 

at a statistical 95% CL. In the DFSZ model this implies nia < 5 meV/ cos^ (3. 
This is the most restrictive limit on the axion-electron interaction. 

6 Supernova 1987A 
6.1 Energy-Loss Argument 

About two dozen neutrinos from SN 1987A were observed almost twenty years 
ago in several detectors [63]. The total number of events, their energies, and 
the distribution over several seconds correspond reasonably well to theoretical 
expectations. In the standard picture [64,65], the core collapse of a massive 
star leads to a proto neutron star, a solar-mass object at nuclear density and 
a temperature of several 10 MeV, where even neutrinos are trapped. The long 
time scale of emission is explained by diffusive neutrino energy transport. The 
emission of more weakly interacting particles can be a more efficient energy- 
loss channel, resulting in a reduced neutrino burst duration. The late-time 
signal is most sensitive to such losses because the early neutrino emission is 
powered by accretion and thus not very sensitive to volume losses. 

This argument has been applied to many cases, from right-handed neu- 
trinos to Kaluza-Klein gravitons, but axions are the earliest and most widely 
discussed example [23,66-75]. They are emitted by nucleon bremsstrahlung 
N + N N + N + a that depends on the axion-nucleon Yukawa coupling 
QaNN, here taken to be an average of the couplings to neutrons and protons. 
Figure [2] illustrates that axion emission leaves the signal duration unchanged 
when QaNN is very small. For larger couplings, the signal is shortened until it 
reaches a minimum, roughly when the axion mfp corresponds to the geomet- 
ric size of the SN core. For even larger couplings, axions are trapped and are 
emitted from an "axion sphere." When it moves beyond the neutrino sphere, 
the signal duration once more remains unaffected. 

Of course, such "strongly" interacting axions are not necessarily harmless. 
They may play an important role during the infall phase. Moreover, in the 

^ The limit on gaee stated in [61] is an order of magnitude more restrictive, but 
this is an obvious misprint. Likewise, their stated hmit rua < 5 meV cos^ f) has 
an incorrect scaling with cos^ (3. 
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Fig. 2. Relative duration of a SN neutrino burst as a function of the axion- nucleon 
coupling [10]. Freely streaming axions are emitted from the entire core volume, 
trapped ones from an "axion sphere." The sohd hne is from numerical calcula- 
tions [71,72]. The dotted line is an arbitrary continuation to guide the eye 



water Cherenkov detectors that registered the SN 1987A neutrinos, these ax- 
ions would have interacted with oxygen nuclei, leading to the release of 7 rays 
and causing too many events [76]. 

However, for axions and other particles, the trapping regime is usually 
excluded by other arguments so that the free-streaming regime is of greater 
interest. An approximate analytic constraint on the energy-loss rate is [9] 

£a < 1 X 10^9 erg g-^s-i , (21) 

to be calculated at g ^ 3 x IQi'^gcm-^ j, ^ 30MeV. If we take the SN 
core to have a mass of about 1 Mq = 2 x 10'^^ g, this corresponds to an axion 
luminosity La = BaM^ = 2 x lO^^ergs"^. The gravitational binding energy 
of the neutron star is about 3 x 10'^'^ erg and the emission lasts up to 10 s, i.e., 
axion losses would compete significantly with neutrino emission. 

We stress that the criterion pT|) is not arbitrary, but was distilled from 
several numerical simulations that consistently showed that the burst duration 
was roughly halved when the limit (PT|) was saturated [9]. Axion losses are 
then not small so that T = 30 MeV is not the unperturbed temperature of 
these models. Different numerical models with different input physics probably 
have internal temperatures that are more similar once significant axion losses 
are included. In any event, (|2ip represents quite accurately the results from 
different simulations. 

Recently, self-consistent cooling calculations were performed for Kaluza- 
Klein gravitons [77] , once more confirming (j2ip . The neutrino signal duration 
was directly compared with the data and the limit, corresponding to pT|) . was 
found to have a 95% statistical CL. 
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Of course, this and any other numerical study rehes on input physics for 
which systematic uncertainties are difBcult to quantify, notably the nuclear 
equation of state and the neutrino opacities. In addition, the data are very 
sparse so that any conclusion based on them suffers from the usual problems 
of small-number statistics. Therefore, ((2T|) should be viewed as a reasonable 
guide as to where a new energy-loss channel causes a significant tension with 
the SN 1987A pulse duration. 



6.2 Axion Emission from a Nuclear Medium 

In order to apply ipij) to axions one needs the emission rate from a hot medium 
at nuclear density. The main emission process is nucleon bremsstrahlung, 
N + N^N + N + a, but a reliable calculation of the rate is difficult. Axions 
couple to the nucleon spin so that bremsstrahlung requires the spin to "jiggle" 
in a collision, i.e., spin-conserving interactions do not contribute. This leaves 
the nuclear tensor force that is only crudely modeled by one-pion exchange 
(OPE). In a dense medium, other problems include the modification of particle 
masses and couplings as well as many-body and multiple-scattering effects. 

One approach to estimate the emission rate relies on linear response the- 
ory where emission, absorption, and scattering of neutrinos, axions and other 
particles depend only on a few "form factors" of the medium, i.e., the dynam- 
ical structure functions [73,78-86]. This approach is perturbative to lowest 
order in the weak interaction between neutrinos (or axions) and nucleons, 
whereas the interactions among the medium constituents are lumped into the 
structure functions. 

Assuming the medium to consist of one species of nonrelativistic nucleons, 
the relevant quantity is the dynamical spin-density structure function [73] 

4 f+°° 

S,{Lu,k)^— dte-*(^(t,k).^(0,-k)) , (22) 

J-oo 

where is the nucleon (baryon) density and cr(t^ k) the spatial Fourier trans- 
form of the nucleon spin-density operator. The basic principles of quantum 
mechanics imply the detailed-balancing condition 

S4-uj,k) = Sa{L^,k)e-'^/^ . (23) 

The structure function obeys the sum rule 

Nb 



_ ^ SAu^, k) = 1 + — ( s, . s, cos(k . r,-) 



(24) 



where s^ is the spin operator of nucleon i. The f-sum rule includes a factor lu 
under the integral and establishes a relation to the average nucleon-nucleon 
spin interaction energy [82]. It is often assumed that also the higher sums 
/ dujuj'^ S{iu, k) exist for all n. 
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The axion absorption rate and the volume energy loss rate are given in 
terms of the structure function as 



where fc = |k| w w is the modulus of the axion momentum. Neutrino scatter- 
ing, emission and absorption rates based on the axial vector current are given 
by similar phase-space integrals. 

A reliable expression for Saii^jk) is not available so that we need to 
use heuristic reasoning. The large nucleon mass compared with the emit- 
ted axion energy suggests to use the long- wavelength approximation 5cr(w) — 
limk^o 5'(t(w, k), i.e., we neglect the momentum transfer to the medium. In 
this limit ([22]) represents essentially the Fourier transform of the autocorrela- 
tion function of a single nucleon spin. 

If we picture the nucleon spin as a classical vector that is kicked by a 
random force, we find [79] 

S^J^^^i.) = , (26) 

where -To- is the spin fluctuation rate. Being a classical result, the quantum- 
mechanical detailed-balancing property is missing. Overall we thus write 

= s{./T) X ( ' tl " ! (27) 



ifTi"^"^^ \e-/^for Lu<0, 

where s{x) is an even function normalized to s(0) = 1. The axion emission 
rate per unit mass, = Qa/e, therefore is 



where 



r°° x^e'"^ r IT 

"-i j^titW*'' '''' 

For Fa/T <C 1 (dilute medium) and assuming s{x) — 1, we find F = FajlT . 

A perturbative calculation, relevant for a dilute medium, and using the 
OPE approximation, yields [70, 79] 

rOPE = A-K^l^o? I^mi = 450 MeV ^ ( ^ ] (30) 

" ^^/2 3xl0i4gcm-3 VaOMeVy ' ^ ' 

where a^^ = (/2mAr/m7r)^/47r w 15 with / « 1.0. For soft energies, brems- 
strahlung depends only on the on-shell spin-dependent nucleon scattering rate. 
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Based on measured nuclear phase shifts it was argued that the OPE resuh 
was an overestimation by about a factor of 4 [75] . 

Either way, T^jT is not small compared to unity, but also not very large. 
A possible range 1 < Fc/T < 10 appears generous. With s{x) = 1 this would 
imply F k, 0.5 for F„/T = 1, a maximum ol F k, 1.35 near Fc /T = 7 and 
F K. 1.3 for Fa/T = 10. Of course, s{x) probably decreases with x or else the 
f-sum and higher sums of Saii-o) diverge. On the basis of existing information 
one cannot do better than assume to be a factor of order unity. 

The SN 1987A limits are particularly interesting for hadronic axions where 
the bounds on aaee are moot. Therefore, we use Cp = —0.4 and C„ = 0. 
Initially the proton fraction is relatively large so that we use Yp = 0.3 to scale 
the emission rate to the proton density. With F = 1 and T = 30 MeV we then 
find £a = 1.4 X 10^^ ergg~^ s~^ so that ([21]) implies 

/a > 4 X 10** GeV and m„ < 16 meV . (31) 

Despite a lot of effort that has gone into understanding the axion emission 
rate, these limits remain fairly rough estimates. 

7 Conclusions 

Astrophysics and cosmology provide the most restrictive limits on the axion 
hypothesis as summarized in Fig. [3l Beginning with cosmology, a cold axion 
population would emerge in the early universe that can make up the dark 
matter, but the required axion mass involves many uncertainties [87]. Galactic 
dark matter axions will be searched by the ADMX experiment in the mass 
range rria = 1-100 )ieV [38,88-90]. 

In addition, a population of hot axions is produced. Before confinement, 
the relevant processes involve quarks and gluons [91,92]. Later the most 
generic process is7r + 7r<^7r + a [93]. Axions decouple after the QCD epoch if 
/a < 3 X 10'^ GeV (ma > 0.2 eV). Some of these hot dark matter axions would 
be trapped in galaxies and galaxy clusters. An unsuccesful search for a decay 
line [94-96] provides direct limits on a range of axion masses marked "Tele- 
scope" in Fig. [3l Moreover, the usual structure-formation arguments provide 
the hot-dark matter limits [97,98]. Axions decay on a cosmic time scale for 
nia ^ 20 eV. The decay photons would cause a variety of observable conse- 
quences [99], seamlessly connecting with the hot dark matter limit so that 
cosmology alone rules out axions in the entire mass range > 1 eV. 

Figure [3] also shows the stellar-evolution limits discussed in this chapter, 
notably the globular-cluster limit on the axion-photon coupling. The globular- 
cluster and white-dwarf limits for DFSZ axions with cos^ /5 = 1 are shown as a 
light-grey exclusion bar. Even in the DFSZ model, the axion-electron coupling 
could be accidentally small and at tree level it is entirely absent for hadronic 
axions. Therefore, these limits are far less generic than those based on the 
axion-photon or axion-nuclcon interaction. 
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The requirement that the neutrino signal of SN 1987A was not excessively 
shortened by axion losses pushes the hmits down to rua ^ 10 meV. However, 
this limit involves many uncertainties that are difficult to quantify so that 
it is somewhat schematic. The CAST search for solar axions [46] covers new 
territory in the parameter plane of nia and ga-y-y, but a signal would represent 
a conflict with the SN 1987A limit. While this limit certainly suggests that 
axions more plausibly have masses relevant for cold dark matter, a single 
argument, measurement or observation is never conclusive. 

In the DFSZ model, the limits from white-dwarf cooling based on the 
axion-electron interaction and those from SN 1987A from the axion-nucleon 
interaction are quite similar. Therefore, axion emission could still play an 
important role as an energy-loss channel of both SNe and white dwarfs and 
for other evolved stars, e.g. asymptotic giant stars. 

In summary, axions provide a show-case example for the fascinating inter- 
play between astrophysics, cosmology and particle physics to solve some of 
the deepest mysteries at the interface between inner space and outer space. 
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